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Abstract 



In this paper, we study the discrete cubic nonhnear Schrodinger 
lattice under Hamiltonian perturbations. First we develop a complete 
isospcctral theory relevant to the hypcrbohc structures of the lattice 
without perturbations. In particular, Backlund-Darboux transforma- 
tions are utilized to generate heteroclinic orbits and Melnikov vectors. 
Then wc give coordinate-expressions for persistent invariant manifolds 
and Fenichel fibers for the perturbed lattice. Finally based upon the 
above machinery, existence of codimension 2 transversal homoclinic 
tubes is established through a Melnikov type calculation and an im- 
plicit function argument. We also discuss symbolic dynamics of invari- 
ant tubes each of which consists of a doubly infinite sequence of curve 
segments when the lattice is four dimensional. Structures inside the 
asymptotic manifolds of the transversal homoclinic tubes are studied, 
special orbits, in particular homoclinic orbits and heteroclinic orbits 
when the lattice is four dimensional, are studied. 



Keywords: Homoclinic tubes, Backlund-Darboux transformations, 
cubic nonlinear Schrodinger lattice, Melnikov vectors, Fenichel fibers. 
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1 Introduction 



The concept of a homoclinic tube was introduced by Silnikov [16] in a study 
on the structure of the neighborhood of a homochnic tube asymptotic to an 
invariant torus a under a diffeomorphism F in a finite dimensional phase 
space. The asymptotic torus is of saddle type. The homochnic tube consists 
of a doubly infinite sequence of tori {cjj , j = 0, ±1, ±2, • • •} in the transversal 
intersection of the stable and unstable manifolds of cj, such that (Tj+i = Fodj 
for any j. It is a generalization of the concept of a transversal homoclinic 
orbit when the points are replaced by tori. Silnikov obtained a similar 
theorem on the symbolic dynamics structures in the neighborhood of the 
homoclinic tube as Smale's theorem for a transversal homoclinic orbit [15] 
[14]. 

We are interested in homoclinic tubes for several reasons: 1. Especially in 
high dimensions, dynamics inside each invariant tubes in the neighborhoods 
of homoclinic tubes are often chaotic too. We call such chaotic dynam- 
ics ^^chaos in the smaU\ and the symbolic dynamics of the invariant tubes 
'^chaos in the large" . Such cascade structures are more important than the 
structures in a neighborhood of a homoclinic orbit, when high or infinite 
dimensional dynamical systems are studied. 2. Symbolic dynamics struc- 
tures in the neighborhoods of homoclinic tubes are more observable than 
in the neighborhoods of homoclinic orbits in numerical and physical experi- 
ments due to the large dimensionality and the robustness of the homoclinic 
tubes. 3. When studying high or infinite dimensional Hamiltonian system 
(for example, the discrete NLS or NLS equations under Hamiltonian pertur- 
bations), each invariant tube contains both KAM tori and stochastic layers 
(chaos in the small). Thus, not only dynamics inside each stochastic layer 
is chaotic, all these stochastic layers also move chaotically under Poincare 
maps. 

In [7], we proved the existence of transversal homoclinic tubes which 
are asymptotic to locally invariant center manifolds for the cubic nonlinear 
Schrodinger equaton under Hamiltonian perturbations. Due to the locally 
invariant nature of the center manifolds, the symbolic dynamics in the neigh- 
borhood of the homoclinic tube is very difficult to establish. The locally 
invariant nature is very difficult to control in infinite dimensions. There- 
fore, we need to seek low dimensionality and hopefully can establish the 
above mentioned symbolic dynamics for the low dimensional systems which 
is the first step toward establishing such symbolic dynamics for PDEs. Then 
naturally we want to study the finite difference discretization of the cubic 
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nonlinear Schrodinger equation under Hamiltonian perturbations (perturbed 
discrete NLS). The dimension of the perturbed discrete NLS is n. As shown 
in this paper, when n = 4, the center manifold is actually invariant and 
Silnikov's theorem implies the symbolic dynamics of line segments. 

Denote by VF^'^^ a normally hyperbolic locally invariant center manifold, 
by and W^'^^^ the locally invariant center-unstable and center-stable 

manifolds such that W^'^^ C W^'^'"^ nW^'^^\ and by F* the evolution operator 
of the system. We call Tl a transversal homoclinic tube if W C W^'^'^^nW^'^'^\ 
the intersection between VF^'^") and l^^*^*) is transversal at H, and H has 
the same dimension with W^^\ Let S be an appropriate Poincare section, 
and P is the Poincare map induced by the flow F*; then H PI S is called a 
transversal homoclinic tube under the Poincare map P. 

In [10], the discrete cubic nonlinear Schrodinger equation under dissipa- 
tive perturbations is studied. Existence of a symmetric pair of homoclinic 
orbits is established through a Melnikov calculation and a geometric argu- 
ment. In [11], symbolic dynamics in the neighborhood of the symmetric pair 
of homoclinic orbits is established. In contrast to these previous works, the 
current work is a study on Hamiltonian perturbations and homoclinic tubes. 

The paper is organized as follows: Section 2 is on the formulation of the 
problem, section 3 is on the isospectral theory for the discrete NLS equation, 
section 4 is on the coordinatization for invariant submanifolds, section 5 is 
on the existence of transversal homoclinic tubes, section 6 is on the symbolic 
dynamics of segments for the case when the perturbed discrete NLS equation 
is 4-dimensional, section 7 is on structures inside the asymptotic manifolds 
of the transversal homoclinic tubes, and section 8 is the conclusion. 
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2 Formulation of the Problem 



Consider the discretized cubic nonlinear Schrodinger equation under Hamil- 
tonian perturbations (perturbed discrete NLS), 

iQn = PndH/dqn , (2.1) 



where 

H = Hq + eHi , 

n=0 

£ is the perturbation parameter, eHi is the Hamiltonian perturbation, i = 
\/— 1 is the imaginary unit, a; is a positive parameter, and p„ = 1 + /i^|gf„p. 
g„ satisfies the periodic and even boundary conditions, 

Qn+N = Qn, Q-n = Qn , (2-2) 

where AT is a positive integer N > 3. (2.1) is a 2(M + l)-dimensional system, 
where M = N/2 {N even) and M = (AT - l)/2 (AT odd). 

Qniln+i + Qn-i)^ itself is also a constant 

of motion. This invariant, together with Hq, implies that Ylin=o^^Pn ^ 
constant of motion too. Therefore, 

N-l 

D^=l[pn (2.3) 

n=0 

is a constant of motion. 
The phase space is defined as 



Remark 2.1 When e = 0, E^"^ 



s = 



-q, q = iqo,qi, ■■■,qN-i) , 



qn+N = qn, qN-1 



= 9n} ■ (2.4) 



In S (viewed as a vector space over the real numbers), we define the inner 
product, for any two points q^^^ and q^'^\ as follows: 

n=0 
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And the norm of q is defined as \\q \\'^ = {q,q)- We also use tlie notation 
Qn = {qn,rnf, where r„ = 

We will study transversal homoclinic tubes in this phase space S. When 
e = 0, the unperturbed Hamiltonian system is the integrable discrete cubic 
nonlinear Schrodinger equation, 

iQn = -^[Qn+l - 2gn + Q'n-l] + knPC^n+l + Qn-l) - ^U^Qn ■ 

The Hamiltonian perturbation term Hi can be very general. In this paper, 
for concreteness, we will study the simple example. 



where a\ and a2 are real parameters. The corresponding equation is. 



We will show that there exist transversal homoclinic tubes for this system 
in the phase space S through Melnikov analysis and an implicit function 
argument. 

Remark 2.2 In the literature [2], for integrable systems under Hamiltonian 
perturbations, an interesting fact is that the Melnikov function sometimes is 
identically zero. Consider the NLS equation under Hamiltonian perturbation 




n=0 




(2.5) 



[2J 



-2\q\^q-£q, 



[XXXX ) 



(2.6) 



with the Hamiltonian 




In [2], the Melnikov function M is built with Hq 



M = 
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where Q is a heteroclinic orbit generated through Bdcklund-Darhoux transfor- 
mations, which is asymptotic to a periodic orbit independent of x in \t\ — oo 
limit. The reason why M is identically zero is as follows 



M 



1 r°° dH 
e J-oo dt 



1 

dt = — 
Q ^ 



H{ lim Q) - H{ lim Q) 



. 



From the study [9], one doubts that Hq is the proper invariant to build the 
Melnikov function, since the center-unstable manifold associated with the 
unstable mode is really a level set of the invariant Fi (3.6). Nevertheless, 
direct calculation shows that 

(Q) = linear combination of { (Q) , (Q) , ^ (Q) , ^ (Q) 1 
oq L oq oq oq oq ) 

with coefficients only dependent upon t , 

where I = Jq \q\'^dx is still invariant under the perturbed flow (2.6). This 
relation and the above Melnikov function do not imply that the Melnikov 
function built with Fi is identically zero, but offer a hint that it is prob- 
ably identically zero too. Prom some trial calculations for this paper, if a 
perturbation term still keeps 

N-l 
n=0 

invariant, the corresponding term of the Melnikov function built from Fi 
(3.6) is identically zero. Such perturbation terms are for example, 

^ N-l 
n=0 

When the Melnikov function is identically zero, second order Melnikov func- 
tion is needed to measure the splitting between center-unstable and center- 
stable manifolds [8j. 
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3 Isospectral Theory for Discrete NLS Equation 

In this section, we will study the isospectral theory for the discrete nonlinear 
Schrodinger equation 

iqn = -^[qn+l - 2gn + Qn-l] + \qnfiqn+l + Qn-l) - ^LO^Qn , (3.1) 

where i = QnS are complex variables, n G Z; under periodic and even 

boundary conditions, 

Qn+N = Qn ) Q-n = Qn ■ 

Since it is brand new, we will present the isospectral theory for the discrete 
NLS equation in details. 

The integrability of the NLS equation (3.1) is proven with the use of the 
discretized Lax pair [1] : 

ifn+l = (3.2) 
ipn = (3.3) 



where 
Li-) = 



z ihqn 
ihqn l/z 



^{z) — j_( I - + 2iXh - h^qnQn-i + uP'h?' -izhqn + {l/z)ihqn-i 

~ \ -izhqn-i + {l/z)ihqn 1/ z^ - I + 2iXh + h'^qnqn-i - oj'^h'^ 

and where z = exp{iXh). Compatibility of the over determined system 
(3.2,3.3) gives the "Lax representation" 

of the discrete NLS equation (3.1). Focusing attention upon the discrete 
spatial flow (3.2), we let Yi^\ Y^"^^ be the fundamental solutions of the ODE 
(3.2), i.e. solutions with the initial conditions: 

The Floquet discriminant: 

A:CxS^C (3.4) 



9 



is defined by 

/^{z-q)=tr{M{N-z-q)} , (3.5) 

where 5 is the phase space defined in (2.4), and M(n; z;q) = columns{yi^'*, yi^^} 
is the fundamental matrix of (3.2). 

Remark 3.1 A.{z;q) is a constant of motion for the integrable system (3.1) 
for any z E C [13]. Since A(z; q) is a meromorphic function in z of degree 

(+N, —N), the Floquet discriminant A{z;q) acts as a generating function 
for {M + 1) functionally independent constants of motion, and is the key to 
the complete integrahility of the system (3.1), where M = N/2 (N even), 
M = {N-l)/2 (N odd). 

The Floquet theory here is not standard as can be seen from the Wronskian 
relation: 

Wjv(V'+,V'~) = ^' Wo(^+,^-), 
where D is defined in (2.3), 

-0^ and ■ip~ are any two solutions to the linear system (3.2). In fact, 
Wn+i{il^~^ = PnWn{ip~^ Due to this nonstandardness, modifica- 
tions of the usual definitions of spectral quantities [6] are required. The 
Floquet spectrum is defined as the closure of the complex A for which there 
exists a bounded cigenfunction to the ODE (3.2). In terms of the Floquet 
discriminant A, this is given by 

(t(L) = |zgC -2D < A{z;q) <2D^ . 

Periodic and antiperiodic points z^ are defined by 

A{z';q) =±2D. 
A critical point z'^ is defined by the condition 

= 0. 



dA 

dz 



A multiple point is a critical point which is also a periodic or antiperiodic 
point. The algebraic multiplicity of is defined as the order of the zero 
of A{z) ± 2D. Usually it is 2, but it can exceed 2; when it does equal 
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2, we call the multiple point a double point, and denote it by z'^. The 
geometric multiplicity of z"* is defined as the dimension of the periodic (or 
antiperiodic) eigenspace of (3.2) at z™, and is either 1 or 2. 
The normalized Floquet discriminant A is defined as 

A = A/D . 

Remark 3.2 (Continuum Limit) In the continuum limit (i.e. h ^ 0), 

the Hamiltonian has a limit in the manner: H/h He, where is the 
Hamiltonian for NLS PDE, He = i jQ{qxQx + 2aj^|gp — The Lax 

pair (3.2;3.3) also tends to the corresponding Lax pair for NLS PDE with 
spectral parameter X (z = e^^^) [9]. IfQ = max„{|g„|} is finite, then pn ^ ^ 
as h ^ 0. Therefore, = {11^=^0^ Pn} 1 as h ^ 0. The nonstandard 
Floquet theory for the spatial part of the Lax pair (3.2) becomes the standard 
Floquet theory in the continuum limit. 

3.1 Examples of Floquet Spectra 

Consider the uniform solution: qn = qc, Vn 

qc{t) = a exp | - i[2{a^ - u^)t - 7] | . 

The corresponding Bloch functions of the Lax pair are given by: 

/ (i - exp{-i[(a2 - u^)t - 7/2]} 



V'n = (Vpe-^^)'^e^-* 
where 



— ■j/iaexp{i[(a^ — a;^)t — 7/2]} 
—ihaexp{—i[{a'^ — ui'^)t — j/2]} ^ 

{z - ^e-^f^) exp{i[(a2 - u^)t - 7/2]} j 



2; = Y^cos /? + Y pcos^ /? — 1 , p = l + h^a'^, 



The Floquet discriminant is given by: 

A = 2Dcos(iV/3) , 
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Figure 3.1: An illustration of the Floquet spectra for the uniform potential 
Qc (The number and locations of the spectral points can vary). 



where D = p^l"^. Thus the Floquet spectra are given by: 

-1 < cos(iV/3) < 1 . 



See figure 3.1 for an illustration of the Floquet spectra. Periodic and an- 
tiperiodic points are given by: 



where Zm is a periodic point when m is even, and Zm is an antiperiodic 
point when m is odd. See figure 3.1 for an illustration of the periodic and 
antiperiodic points. The following facts are obvious, 

• If A'^ is odd, all the periodic and antiperiodic points are on the real 
axis for sufficiently large |gc|- 

• If A^ is even, except the two points z = ±i, all other periodic and 
antiperiodic points are on the real axis for sufficiently large \q^. 

Derivatives of the Floquet discriminant A with respect to z are given by: 




dA/dz = 2ArL>sin(A^^)[z^sin/3]-ypcos2^- 1 , 
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N cos (AT/?) sin (3 {p cos^ (3-1) 



+ (1 -p) cosj3si-a{Nf3) + y/psm{Nj3) srn^ [3 ^ p cos^ f3 - I . 

The critical points are given by: 

/3= ^vr (/3^0,7r), or cos^ /? = 1/p . 

See figure 3.1 for an illustration of the critical points and double points. 
There can be multiple points with algebraic multiplicity greater than 2. For 
example, when two symmetric double points on the circle collide at the 
intersection points z = ±1, we have multiple points of algebraic multiplicity 
4. In such case, pcos^ P — 1 = and sin(Ar/3) = 0; then (PA/dz"^ = 0, at 
z = ±1. 



3.2 An Important Sequence of Invariants 
Definition 1 The sequence of invariants Fj is defined as: 

Fj{q)=A{zf{q);q). (3.6) 

These invariants Fj^s are perfect candidate for building Melnikov functions. 
The Melnikov vectors are given by the gradients of these invariants. 

Lemma 3.1 Let Zj^\q) be a simple critical point; then 



p(.1 = |^(#^(.l;.r (3.7) 

Oqn OQn ■' 



n+l 



(3.8) 



where ip^ = {4'n^'^\'4'n^'^^)'^ o-f^ two Bloch functions of the Lax pair (3.2,3.3), 
such that 

where t/^^ are periodic in n with period N, Wn = det {tp^,ilj~). 
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Proof: By the definition of critical points, 

A'(zf^(g);g1 = 0. 



Differentiating this equation, we have 



5z 



1 6A' 



Since Zj'^\q) is a simple critical point, z^'^^ is a differentiable function. Thus 



6F^ 



6A 

5A_ 

Sqn 



dA 

,(c) ^ 'dz 



ic) 



z("> Sqn 



This proves equation (3.7). Next we derive the formula (3.8). Let M„ be 
the fundamental matrix to the system (3.2), i.e. the matrix solution to (3.2) 
with initial condition Mq being a 2 x 2 identity matrix. Variation of q leads 
to the variational equation for the variation of M„ at fixed z, 

^ ih Sqn ^ 

SMn + 



SM, 



n+l 



5Mn = . 




ih 6qn 







Mr, 



Let 6Mn = MnAji, where A„ is a 2 x 2 matrix to be determined, we have 

(3.9) 



where 



SUn 




Solving the system (3.9), we have 
SMn = Mn 
SMo = . 
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Then, 



Thus, 



5A{z,q) = trace |Mjv[Xj M,"^5C/j_iM,_i]| . 

^ = ih trace |m-_^\ ( ) ^"^^} ' ^^"^O) 
^ = -ih trace |m-^, ( 1 ] ^"^^} ' (^-^l) 



where rn = —q,i- Let -0"^ and ip be two Bloch functions for the discrete 
Lax pairs (3.2,3.3), 

i^t = I)"/^C^"/^^^ , (3.12) 

where ip^ are periodic in n with period A'^. Let be the 2x2 matrix with 
and as the column vectors, 

5„ = (V^+,V'-) . (3.13) 

Then 

M„ = S^Bo"^ , Mn = Bo(^^^ ^^„? ^ Bo~^ . (3.14) 
Substitute the representations (3.13) and (3.14) into (3.10,3.11), we have 



Sqn 5qn 2Wn+l I (+,i) 



where 

Wn = det Bn, Vn = ■ 

Thus, 



.5A _ ihjC - C') 

Sqn 2Wn+l 



which is the formula (3.8). The lemma is proved. □ 
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3.3 Backlund-Darboux Transformations 



The hyperbolic structure and homochnic orbits for (3.1) are constructed 
through the Backlund-Darboux transformations, which were built in [6]. 
First, we present the Backlund-Darboux transformations. Then, we show 
how to construct homoclinic orbits. 

Fix a solution Qnit) of the system (3.1), for which the linear operator 
Ln has a double point z"^ of geometric multiplicity 2, which is not on the 
unit circle. We denote two linearly independent solutions (Bloch functions) 
of the discrete Lax pair (3.2;3.3) a.t z = z'^ by {(j)^,4)~). Thus, a general 
solution of the discrete Lax pair (3.2;3.3) at {qn{t), z"^) is given by 

(t)n{t- z'^, c+, c~) = c+<;^+ + c~(i)~, 

where c"*" and c~ are complex parameters. We use ^„ to define a transfor- 
mation matrix r„ by 



where. 



an = 



dn = 



Z + {l/z)an bn 

Cn -l/z + zdn 



M + \z''\^\<l>nl? 



1 



Z'^Ar 



M + \z''\^\4>nl 



1 



A„, = 



{Z'^f^n 

\z'^\'^-l 

Z'^z'^t^n ' 
1 



2 I ud|2| 



\(l)nir + \z' 



!>n2 



From these formulae, we see that 

ciji = dji, hji = Cyi- 
Then we define Qn and by 



(3.15) 
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and 



(3.16) 



where jpn solves the discrete Lax pair (3.2;3.3) at {quit), z). Formulas (3.15) 
and (3.16) are the Backlund-Darboux transformations for the potential and 
eigenf unctions, respectively. We have the following theorem [6]. 

Theorem 3.1 (Backlund-Darboux Transformations) Letqnit) denote 

a solution of the system (3.1), for which the linear operator Ln has a dou- 
ble point z'^ of geometric multiplicity 2, which is not on the unit circle and 
which is associated with an instability. We denote two linearly independent 
solutions (Bloch functions) of the discrete Lax pair (3.2;3.3) at {qn,z'^) by 
(^+,^-). We define Qn{t) and 'Sn{t;z) by (3.15) and (3.16). Then 

1- Quit) is also a solution of the system (3.1). (The eveness of Qn can 
be obtained by choosing the complex Bdcklund parameter c'^/c" to lie 
on a certain curve, as shown in the example below.) 

2. '^n{t'-,z) solves the discrete Lax pair (3.2;3.3) at {Qn{t),z). 

3. A{z- Qn) = A{z; g„), for all zeC. 

4- Qn{t) is homoclinic to qn{t) in the sense that Qn{t) e*^± qn{t), 
exponentially as exp(— cr|t|) as t ^ ±00. Here 9± are the phase shifts, 
a is a nonvanishing growth rate associated to the double point z'^, and 
explicit formulas can be developed for this growth rate and for the phase 
shifts 6±. 



This theorem is quite general, constructing homoclinic solutions from a wide 
class of starting solutions qn{t). Its proof is by direct verification [6]. 



3.4 Heteroclinic Orbits and Melnikov Vectors 



Through Backlund-Darboux transformations, heteroclinic orbits Qn are gen- 
erated, with the formula expression (3.15) 

Qn = ih ^bn+1 — ttn+iqn 

ih-\l ■ 



+ Z'''><ln{\^): 



+11 
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The Melnikov vector field located on this heteroclinic orbit is given by 



I 



SQr. 



-2 ^(1) ^(1) \ 



(3.17) 



where 



w„.= 



En 



1)|2 ^ Ud)\2 



K 



where 



3.4.1 An Important Example 

In this subsection, we start with the uniform solution of (3.1) 

9n = 9c , Vn; gc = aexpj - i[2(a^ - a;^)i - 7]| . (3.18) 
We choose the amplitude a in the range 



AT tan — < a < A/" tan — , when A/" > 3 , 
N N 



TT 

3 tan — < a < oo , when AT = 3 ; 



(3.19) 



so that there is only one set of quadruplets of double points which are not on 
the unit circle, and denote one of them by 2; = = z^^ which corresponds 
to (5 = tt/N (see subsection 3.1 and figure 3.1). The heteroclinic orbit Q„ 
is given by 



Qn — Qc{En+l] 



An+l - 2 COS (3 J p COs2 P - lBn+1 



(3.20) 
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and the Melnikov vectors evaluated on this heterochnic orbit are given by 



SFi 



K 



sech [2nt + 2p] i ^^(2) 



(3.21) 



where 



En = hacos (3 + 1 sech + 2p] cos[(2n - l)/3 + 1?] , 

An+i = hacos P + ^/9Cos2/3- 1 sech [2/xt + 2p] cos[(2n + 3)/3 + , 
Bn+i = cos if + i simp tanh[2/xt + 2p] + sech [2/xt + 2p] cos[2(n + l)/3 + , 



1(1) = 



cos /3 sech [2fit + 2p] + cos[(2n + 1)P + 'd + (p] 



-i tanh[2//i + 2p] sin[(2n + I)f3 + 'd + ip] 



w = 



COS P sech [2/xt + 2p] + cos[(2n + l)/3 + 7? - t/p] 



-ztanh[2/Ltt + 2p\ sin[(2n + l)/3 + - 



i2e(t) 



K = -2Nh'^a{\ - z^)[8p3/2z2]-ypcos2/3-l , 



P = 'k/N, p = l + h^a^ , fi = 2h~^y/psmf3^Jpcos'^f3-l 
h = l/N, c+/c^ = ie^Pe'^ , i? G [0, 27r] , pe (-00, 00) , 
z = y^cos /3 + yjpcos^ p-1 , 6'(t) = (a^ - a;^)i - 7/2 , 



' pcos"^ /? — 1 + i-y/psin (3 = hae^^ , 

where if = sin~-^[y^(/ia)~"^ sin/3], ip G (0, 7r/2). 

Next we study the "evenness" condition: Q_„ = Qn- It turns out that 
the choices ?? = — /? , — /3 + tt in the formula of Qn lead to the evenness of 
Qn in n. In terms of figure eight structure of Qn, "& = —f3 corresponds to 
one ear of the figure eight, and i3 = —f3 + vr corresponds to the other ear. 
The even formula for Qn is given by. 



Qn — Qc 



r/A„ - 1 



(3.22) 
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where 

r = 1 — cos 2(f — i sin 2ip tanh[2/xt + 2p] , 

An = 1 lb cos(/7[cos (3]^^ sech[2^t + 2p] cos[2nf3] , 

where ('+' corresponds to i} = —(3). The Melnikov vectors evaluated on 
these heterochnic orbits are not necessarily even and are in fact not even 
in n. For the purpose of calculating the Melnikov functions, only the even 
parts of the Melnikov vectors are needed, which are given by 



5Fi 



k^"^ sechp^i + 2p][n„]-^ ( ^I* 



even 



-X, 



(2,e) 



where 



(3.23) 



cos (3 ± cos LP sech[2//i + 2p\ cos[2(n — l)f5] 
cos (3 lb cos if sech[2/xt + 2p] cos[2(n + 1)(3\ 



X. 



cos (3 sech [2//i + 2p\ zb (cos ip 
—isnup tanh[2/xt + 2p] ) cos [2n/3] 



J26{t) 



cos (3 sech [2//i + 2p\ zb (cos Lp 
+i sin (p tanh[2/xt + 2p] ) cos [2n/3] 



The heteroclinic orbit (3.22) represents the figure eight structure as il- 
lustrated in Fig.3.2. If we denote by S the circle, we have the topological 
identification: 

(figures) ® 5= IJ Qn{p,l,a,uo,±,N) . 

p£(—oo,oo), 7G[0,27r] 
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Figure 3.2: A geometric illustration of the figure eight structure and its 
corresponding spectral identification. 

4 Coordinate-Expressions for Invariant Submani- 
folds 

In this section, we will give expressions for invariant submanifolds in coor- 
dinates, which enable us to do a Melnikov analysis and an implicit function 
argument. 

4.1 Coordinate-Expressions for Lineeir Invciriant Submani- 
folds 

Consider the discrete cubic integrable nonlinear Schrodinger equation (3.1), 
denote by qc the uniform Stokes solution, 

= oe^^W , e{t) = -[2(a2 - u;^)t - 7] . 

Let 

and linearize equation (3.1) at Qc, we have 

^^n = To i^n+l - 2g„ + qn-l] 
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+ a^[qn+i + qn-i] + 2a^^n ■ 
Assume that takes the form, 



COS kju , 



where kj = 2jn/N, (j = 0, 1, • • •, M). Then, 



fi}^^ = ±2smkj Va^+N^ya"^ - tan'^[kj/2] . 

In this paper, we only study the case that a lies in the range given in (3.19), 
so that only are real and nonzero. In fact, 0,^^^ are zero, and ^j^^ are 
imaginary for {j > 1). When j = 1, 

g„ = aS±)e^i*'*e±*'^cosfein, 
where are real constants; 

■& = -^arctan|[(a^ + iV2)sin2(27r/Ar) 

(a^ - NHan\Tr/N))]^/^[N^ - {N^ + a^) cos{27r/N)]-^ 
Denote by B the block in the phase space S, 



o»7 



a + {he'^ + 626" ) cos kin 
, where a e (Artan[7r/Ar], iVtan[27r/iV]), 7 G [0,27r); 



M 

+ ^ Cj COS kjU 

i=2 

61,62 are real; Cj is complex; and i? is given above| . 

In terms of the coordinates {a, 7, 61, 62, Cj(2 < J < M)}, the linear invariant 
center manifold C^'^^ is given by. 



61 = 62 = L 



the linear invariant center-unstable manifold C^'^'^^ is given by. 



the linear invariant center-stable manifold C^^^^ is given by. 
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4.2 Coordinate-Expressions for Persistent Locally Invariant 
Submanifolds 

Under the perturbed flow (2.5), the hnear invariant submanifolds C^^\ 
and C^'^^^ perturb into locally invariant submanifolds W^'^\ PF^*^") and W^'^^ 
(local invariance means that orbits can only enter or leave the submanifolds 
through their boundaries). For references on the proof of such results, see 
for example [4] [10] [12]. In terms of the coordinates {a, 7, 61, 62, Cj (2 < j < 
M)}, for any small (5i > 0, there exist 5i > {I = 0,2,3), such that W^"^ 
has the expression, 

\ bi = V) 7, c; e, a, w) , 

l^(cw) tiie expression, 

h2 = /^"^(a,7,6i,c;e,Q;,a;) ; 
and W^'^'^^ has the expression, 

^'i = /^*^(a,7,62,c;£,a,a;) ; 

where //"^ (/ = 1,2), /(") and /(^) are C" smooth functions for some large 

n, a = (ai,Q!2), |e| < 5o, 7 G [0,27r), \bk\ < ^2, {k = 1,2), lo G Es,, 
a £ Es^, Es^ = (iVtan(7r/iV) + (^1, tan(27r/iV) - Si) when > 4 and 
= (A^tan(-7r/A^) + Si, A) when = 3, 4, where A is a fixed large constant, 
c = (c2, • • -jCm), \\cj\\ < S3 (2 < j < M). Denote by A the annulus, 

A=^qeB 61 = 62 = c = 0, a G 7 G [0,27r)| . (4.1) 

Then, A C >C(^), and A C W^"^; thus, 

//')(a,7,0;£,a,a;) =0, (Z = l,2). 
Notice also that W^"^ = M^^^^") n then. 
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4.3 Coordinate-Expressions for Fenichel Fibers 

The center-unstable and center-stable manifolds W^'^'^^ and W^'^^^ admit fi- 
bration through Fenichel fibers. For references on the proof of such results, 
see for example [10] [12]. Take unstable Fenichel fibers, for example, they are 
a family of curves in PT^'^"), each of them is labeled by a base point in W^'^\ 
is a union of the family of curves over W^^l In this sense, Fenichel 
fibers are coordinates for W^'^'^\ Indeed, they are very good coordinates. 
Fenichel fibers depend smoothly on their base points. In backward 

time, orbits starting from points on the same fiber (in particular, from the 
base point) approach each other exponentially. 

The unstable Fenichel fibers {J^q^^} have the coordinate-expressions: 

a = a^")(6i;ao,7o,Q);£,Q;,w), 

7 = 7(")(6i;ao,7o,co;e,a,a;), 

= c (")(6i;ao,7o,co;e, a,a;), 

b2 = /W(aW,7^"\&i,c(");£,a,a;); 

which are smooth in bi, and C"~^ smooth in (ao,7o,Q)) and {e,a,u)), 
where hi G {—62,62). The base points q are given by, 

a(^\f[''^;ao,jo,oo;£,a,u) = ao , 

7^"-'(/i'''';ao,7o,co;e,a,u;) = 70 , 
c ^"H/i''';ao,7o,co;e,a,a;) = cq , 
/H(aH,^{«),/W,c-H;e,a,u;) 

= /^"-'(ao,7o,/f'',co;e,a,u;) 

= /2''''(«o,7o,co;e,a,a;) , 

where /[""^ = (oq , 70 , cq ; e, a, tj) . 

The stable Fenichel fibers {J^q^^} have the coordinate-expressions: 

a = a^^\b2;ao,^o,oo;e,a,uj), 

7 = 7^*''(^'2;ao,7o,co;£,a,a;), 

S = c^^\b2;ao,-fo,co;e,a,uj), (4.2) 

61 = /(^)(a(^),7(^\62,c(^);£,a,a;); 
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which are C" smooth in 62, and C"' ^ smooth in (ao,7o,co) and {s,a,u) 
where 62 € (— ^2,(^2)- The base points q are given by, 



7 



''''\f2^;ao,7o,oo;£,a,u) =70 , 
/(^)(aW,7(^\/i^\cW;£,a,a;) 
= fi'\ao,lo,co;e,a,uj) , 



where f^""^ = f^''\ao,-fo,co;e,a,u}). 



25 



5 Existence of Transversal Homoclinic Tubes 



In this subsection, we will establish the existence of transversal homoclinic 
tubes, based upon the coordinatization in the previous subsection, a Mel- 
nikov function calculation, and an implicit function argument. 
Define the region U as follows. 



h e (^^2,(^2) 

There exists T > such that F'^ oW^^'^^ intersects the region U (where -F"^ is 
the evolution operator of the system (2.5)), the intersection {F"^ oW^'^^^)nU 
has the expression, 

&i = /T"^(«,7,^'2,c;£,a,a;) . 

Since for any fixed finite T, is a C" diffeomorphism, /^"^ is also C" 
smooth. Notice also that when e = 0, W^'^'^'^ = W^'^^\ we have 

Define the function. 



A{ao,'yo,co,b2;e,a,u) = /|"''(a(*\ 7^*), 62, a, w) 

- /W(aW,7(^),fe2,c(^);e,a,a;) , 



where 

a(*) = a^^\b2;ao,jo,oo;e,a,uj), 

^{s) ^ j^^\b2;ao,jo,oo;e,a,uj), 

^(s) ^ c^^\b2;ao,-fo,oo;e,a,uj); 

are defined in the coordinate-expression (4.2) for the stable Fenichel fibers, 
and A is a C"~^ smooth function. Setting 62 = ^62, for sufficiently small e 
and \\co II (||co || = [E,^2 ]'/'), we have 

3 

A(ao,7o,co, -(52;e,a,a;) = sMp^{ao,-fo;a,uj) + eRi{ao,-fo,co;e,a,uj) 

+ e^i?2(ao,7o,Q);£,a,w) , 
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where M^^ (ao, 70; a, w) is the Melnikov function, 

Ri^O{\\co\\) , as ||q)|| ^0 . 

For a derivation on this result, see for example [10]. Define the function A 
as follows. 



= A(ao,7o,Q);£,a,w) 

= -A(ao,7o,co, -c)2;£,a,t^) 

= -^_F^(ao,7o;a,'^) + -Ri(ao, 70, cq; e, a, w) 

+ e-R2(ao,7o,co;e,a,a;) . 



(5.1) 



Then A is C" ^ smooth in (ao, 70, cq; £, a, oj). The Melnikov function Mp^ = 
(a, 7; a, w) is given as: 

Mh=^ {Fi^H,}\^dt = 2 Im {g^Pn^} dt , 

J— c» n=0 I- n 'Jin ) Q 



where 



n=0 



dFi dHi _ dFi dHi 

OQn OTn OTn OQu 



and Q is given in (3.22) and grad Fi{Q) is given in (3.23). Thus, 

N-l 



1 



n=0 
AT-l 



^ In pn + {qn + qn)qn 



dt 



Q 



+ a2 



/oo ^ ' ^ 

-~ n=0 



Im 



dFi 

dr„. 



'^Qn^ lnp„ + {ql + g^)g„ 



d<5.2) 



Setting Mp^ = in (5.2), we have an equation of the form, 

ai — 'koKa2 = , 

where 



(5.3) 



K = /c(a, 7; io, N) 



1 f2{a,r,o^,N) 
4uj /i(a,7;w, AT) 
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Figure 5.1: The graph of n for a fixed value of ui. 



Figure 5.2: The graph of k, when a = oj. 
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Figure 5.3: The graph of |^ for a fixed value of lo. 

The graphs of k arc shown in figures 5.1 and 5.2. Figure 5.2 is for the case 
when Qc is a circle of fixed points. For the convenience of later argument, 
we denote Mp,^ simply by M and denote the surface defined by (5.3) by Sy. 

: 7o = r(o)(ao;a,a;) . (5.4) 

Theorem 5.1 There exist a positive constant Eq > and a region £ for 
(a,Lo), such that for any e G (— eo^eo) o-nd any (a,u;) G £, there exists a 
codimension 2 transversal homoclinic tube asymptotic to the codimension 2 

center manifold W^'^\ 

Proof: There exists a region V on the surface S'-y in (5.4), such that 

d 

-— M(ao,7o;a,^^) 7^ 

and is bounded. See figures (5.3) and (5.4) for the corresponding graphs 
of Next we want to solve the equation (5.1) by the implicit function 

theorem [[3],p265]. For any (oq, 79; a*, a;*) G V, 

^A(a^,7o*,0;0,a*,a;*) = ^M(a^, 70*; a*, a;*) ^ (5.5) 

C70 070 

and is bounded. Then by the implicit function theorem [[3],p265], there is a 
neighborhood W^*) of (ag, Q) = 0; e = 0, a* ,uj*) and a unique C""^ function, 

70 = r(*)(ao,co;£,a,a;) 

defined in W'^*\ such that 

rW(a^,co = 0;£ = 0,a*,a;*)=7o* , 
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Figure 5.4: The graph of |^ when a = lo. 

and 

A(ao,rW(oo,co;e, a,w),co;£,a,w) = . 
Since A is a C"~^ smooth function, by relation (5.5), we have 
d 

—A{ao,T^*\ao,co]e,a,uj),oo;e,a,uj) 7^ , 

(770 

and is bounded for (oq, cq; e, a, cj) G W^*\ Thus, the center-unstable mani- 
fold W^^'^^ and the center-stable manifold W^^^^ have a transversal intersec- 
tion at the neighborhood W^*^ . Let 

(a5,7*;a*,a;*)eV 

then there is a unique C""^ function, 

70 = r(ao,q);£,Q;,w) (5.6) 

defined in W, such that 

r(ao,co = 0;£ = 0,Q;,a;) = T^°\ao;a,oj) , 

and 

A{ao,T{ao,co;e,a,uj),co;e,a,uj) = . 

Notice that (5.6) defines a codimension 1 submanifold W^'^^ of the center 
manifold 1^^"^^ which has the expression, 

bi = fi'\ao,r{ao,co;e,a,u),co;e,a,uj) , 
70 = T{ao,co;e,a,uj) ; 
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Figure 5.5: An illustration of the homoclinic tube. 



where (ao, Q); e, a, w) € W. Define H as follows, 

n= U F*o u .^(^), 

te(-oo,oo) qew^"'^ 

where F* is the evolution operator for the system (2.5). Then H is the 
codimension 2 transversal homoclinic tube asymptotic to W^'^\ □ 

For an illustration of the homoclinic tube, see figure 5.5. Studies on the 
symbolic dynamics in the neighborhood of this homoclinic tube are topics 
of future works. 
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6 Symbolic Dynamics of Segments for the Case 
(A^ = 3, Nonresonant Region) 



n = iqes 



g Vn > . (6-1) 



Denote by 11 the plane, 



n is invariant under the flow governed by (2.5). The annulus A deflned in 
(4.1) is a subset of 11. Away from the resonant circle: 

S^ = i^qeB 6i = 62 = c = 0, a = 7 G [0,27r)| , (6.2) 

the dynamics is given by periodic motions, and denote by A such an annulus 
with boundaries which are periodic orbits. When N = 3, the Melnikov 
analysis shows that there exists a transversal homoclinic tube asymptotic to 
n. By the invariance of the Hamiltonian H, each orbit inside the homoclinic 
tube, approaching a periodic orbit in either the forward or the backward 
time, has to approach the same periodic orbit in both forward and backward 
time. Thus, there exists a homoclinic tube asymptotic to the annulus A, as 
illustrated in figure 6.1. 

Definition 2 Define a Poincare section S as follows, 



qe B 



7 = 



and let P he the Poincare map, 

P : P C S H^- S , 
induced by the flow, where V is the domain of definition for P. 

The transversal homoclinic tube T-Cp for the Poincare map P is illustrated 
in figure 6.2. This transversal homoclinic tube Tip consists of segments, 
and is asymptotic to a segment s in the annulus A. The segment s consists 
of fixed points of P. Invariant tubes in the neighborhood of Hp also consist 
of segments. Wc have the following symbolic dynamics theorem for the 
homoclinic tube Hp. 
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Figure 6.1: A homoclinic tube asymptotic to a nonresonant annulus when 
= 3. 



Figure 6.2: A homoclinic tube asymptotic to a segment under the Poincare 
map P when N = 3. 
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Theorem 6.1 (Silnikov [16]) When N = 3, the set of invariant tubes of 
the Poincare map P lying wholly within a sufficiently small neighborhood of 
the transversal homoclinic tube Hp is in one-to-one correspondence with the 
set of all doubly infinite sequences, 

J = {■ ■ - J-iJoJi,- ■ ■) , 
where ji > for all I and some large integer j* . 
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7 Structures Inside the Asymptotic Manifolds of 
the Transversal Homoclinic Tubes 



Structures inside the transversal homoclinic tubes can be very complicated. 
Since we are studying near integrablc Hamiltonian systems, the asymptotic 
manifold (i.e. the center manifold) VF*-^-* often contains both KAM tori and 
stochastic layers. Studies on such structures are topics of future works. In 
this section, we study interesting structures on an invariant plane inside the 
center manifold W^'^\ and special orbits inside the transversal homoclinic 
tubes. 

The dynamics on 11 (6.1) is governed by 

iq = 2[\q\'^-Lj^]q + el^[aiiq + q) + a2{q'^ + f)]q+[ai+2a2q]^lnp^ , (7.1) 

where p = 1 + h'^\q\'^. Let q = /e*^, where / is the modulus of q; then, 
dl/dt = — esin^ [ai + 4a2-f cos^]^ Inp , 

d^/dt = -2[/2 -w^] -£|2ai/cos^ + 2a2l^cos2^ 



+ 



cos£ 

ai — h 2a2 cos 2^ 



P 1 



The dynamics away from the neighborhood of I = u; is given by periodic 
motions. The dynamics in the neighborhood of J = a; is very interesting. 
Let I = u + rjy, rj = -^/e, and define the resonant annulus A, 



^ G [0, 27r], y is bounded 



Inside the resonant annulus A, the dynamics is governed by, 

dy/dr = fi , 



(7.2) 



/2, 



where 



/i = - sin^ [ai + 4a2(a; + rjy) cos^]^ Inp , 



/2 = -4ujy -rii^2y'^ + 2ai{u} + rjy) cos ^ 

"«! cos ^ 



+2a2 cos 2^ (w + vy) + 



uj + rjy 



+ 2a2 cos 2^ 
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where t = rjt. 



7.1 First Order Dynamics in the Resonant Annulus 

Setting 77 = in (7.2), we have the system, 

dy/dr = — fisin^ [ai + 4a2W cos ^] , 



(7.3) 



d^/dT 



where Q. = h '^p^hipQ and po = 1 + h'^i^'^- This first order system (7.3) is 
also a Hamiltonian system. 



dy/dr = 



dH 



where H = 2a;y^ + Q.{ai cos ^ + 02 w cos 2^). The fixed points of the system 
(7.3) are as follows: 



• When 



4a2W 

i.3) 



< 1, there are four fixed points, 



yy^=0, (j = 1,2,3,4) ; 



(7.4) 



0, ef^=7r, er = -?, 



(3) 



arccos 



When 



> 1, there are two fixed points, 

yS^'^=0, (j = l,2) ; 



The phase diagrams are shown in figure 7.1 when 



"1 



(7.5) 



< 1, and in figure 



7.2 when 



> 1. 
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Figure 7.1: Phase diagrams in the resonant annulus when ^"^^ < 1, (a). 

a2 > 0, ai > 0; (b). 02 > 0, ai < 0; (c). a2 < 0, ai > 0; (d). 02 < 0, 
ai < 0. 



Figure 7.2: Phase diagrams in the resonant annulus when 
ai > 0; (b). ai < 0. 



> 1, (a). 
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7.2 The Full Dynamics in the Resonant Annulus 

First we prove the existence of fixed points for the system (7.2) governing 
the full dynamics in the resonant annulus. 

Lemma 7.1 For any large A > 0, any small Sq > 0, there exists rjQ > 0, 
such that there exist fixed points for the system (7.2), 

y{j) =y0)(^r],ai,a2,uj) , ^^^^ = ^^^\ri,ai,a2,uj) ; 

which are smooth, where (j = 1,2,3,4; when 

> l + So), \r}\ < rjQ, 5o < \aj\ < A, (j = 1, 2), Sq < \u\ < A; 
and 

^^^\0,ai,a2,Lo) = ^i'^ = ^i'\aua2,u;) . 



< 1 — (Jo j and (j 



1, 2; when 



Proof: We want to solve the equations 

/i(y,^;?7,ai,a2,w) = , f2{y,^;v,oii,a2,u) = ; 
in the neighborhoods of the fixed points 

Notice that 

^{y,C;0,ai,a2,uj) = , 
dy 

^{y,C;0,ai,a2,uj) 



Po , 

^Inpo 



4a2a; sin^ ^ — cos ^ (ai + 4Loa2 cos ^) 



^(y,^;0,ai,a2,a;) = -4a; , 

df2 

-^{y,^;0,ai,a2,u}) = . 



(7.6) 
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Then for any 5q > 0, 



det 



4a; po 

/l2 



Inpo 



Aa20J s\T? ^^^^ - cos ^^^^ {ai + 4a2U cos ^^■'^ ) 



when 



4aj«2 



< 1 — 5n or 



4a)a2 



>l + 6o. Then 



dy V^O J 



define Hnear homeomorphisms. Notice that /i and /2 are functions; then 
by the imphcit function theorem [[3], pp265], for any fixed value of Xq"^ 
which satisfies the restrictions given in the Lemma, there is a neighborhood 
Vy'^ of y = {rj = 0,ai,a2,uj), and unique functions 

y^^) = y(^\r],ai,a2,uj), ^^^^ = ^^^)(7?,ai,Q;2,a;); (7.7) 

which are functions defined in Vy^ , such that 

fi{y^^\^^^^;V,ai,a2,u^) = 0, (1 = 1,2) . 

Let V'^^^ = [Jy^y ^'-: then the functions (7.7) are uniquely defined in V'^^'^. 
There exists r/o > 0, such that \rj\ < rjo and the rest of restrictions in the 
Lemma define subregions 
In fact, we have 

and y^^^ {j = 1,2) satisfy the equations. 



7-2 2,2 
I — U + T] 



0, 



where I = u + rjy, '+' for j = 1 and '— ' for j = 2; 



^(3) ^ _^{4) ^ g^j.gp(jg 



ai 
40:2/ 
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and y^^^ = y^^^ satisfies the equation, 



(1 - rj^a2)r - {cj^ + 



2 2 

2 , V af. 



8a2 



??^a2-p- Inp = 



where I = to + rjy. Approximate expressions of these fixed points can be 
obtained: 

y0) = r?yp) + O(r?2) , (i = l,2); (7.8) 



where 



(1) _ 



2/i 

(2) 



1 

4a; 
1 

4uj 



2u{a2UJ + ai) + (2a2 + Q;i/<^)0 Inpo 
2u>{a2UJ — ai) + (2a2 - ai/u))^hipo 



When 



4aJa2 



<1, 



y(0=r,yP + O(r?2) , (/ = 3,4); 



(7.9) 



where 



(3) (4) 1 /Q 2 2 , 2n , «2 /?0 , 



yi-vi 



16a2io 



2oj h? 



^(3) = _^(4) = arccos 



so V 



4a2W 



-1/2 



where if) = . • Let X^-'^ = (y*^-'),^(-'^;?7,Q!i,a2,a;) denote these fixed 



points and notice that X^^ = (yo''\ Co"'^! ^ = 0,01,02,0;) a notation used 
before in (7.6). 

Lemma 7.2 The fixed points X^^ have the same saddle or center nature 



as X, 







Proof: 



dfi 
dy 



sin ^ [4a2 cos $]j^lnp 
sin^ [oi +402(0; + Vey) cos ^]^^^ In p) , 
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sin ^ [4a2 Vs cos (]j^\np 
+^/ssm$, [ai + Aa2{uj + s/ey) cos ^] 



|2(a; + Viy) 



+ 



In p 



Then, 



W + ^/£y 



(7.10) 



in fact, when j = 1,2, ^(X^-^'^) = ^(^^■''^) = 0. Linearize (7.2) at X^^\ 



we have 



where 



1? 



(y\ 




(y\ 















(7.11) 



L3 = 



where 



(/ = 1,2) 



(7.12) 



The eigenvalues of Lj satisfy 



trace(Lj)A + det(Lj) = 0. 



By relations (7.10) and (7.12), 
trace (Lj) = , 

det(Lj) = det 



/f(xf)) f(x(^))\ 



+ 0{'n) ; 



thus, when r] is sufficiently small, the fixed points X^-^) have the same types 
of stabiUty as x'^^\ □ 



41 



The phase diagram is given by the level sets of the rescaled Hamiltonian, 



H = 



Po 



H 



2N 



Po 



where H is the restriction of the Hamiltonian to the resonant annulus A, 



H 



2N ( 



[{Lu + rjyf - 01np] 



+if [ai {cj + riy) cos ^ + 02(0; + r/y)^ cos 2^] In p| . 

H is smooth in rj and has the approximate expression, 

H = 2uy'^ + [ai cos ^ + 020; cos 2^] In po + 0{r]) . 

Setting rj = 0, we have the phase diagrams as shown in figures 7.1 and 
7.2. When r] ^ 0, invariant manifolds of the saddles perturb smoothly, and 
are given by level sets of H; thus, the figure eight loops do not break. By 
Lemmas 7.1 and 7.2, in the resonant annulus A, the phase diagrams when 
rj ^ are topologically equivalent to those when 77 = 0. 



7.3 Special Orbits Inside the Transversal Homoclinic Tubes 

As a consequence of Theorem 5.1, there exist orbits asymptotic to fixed 
points in the resonant annulus A in forward or backward time. See figures 
7.3 and 7.4 for some examples. Whether or not these orbits are actually 
homoclinic orbits or heteroclinic orbits is a very difficult open question. 
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Figure 7.3: Special orbits asymptotic to fixed points in the resonant annulus 



in forward or backward time when 



ai 



4a2W 



< 1, (a). a2 > 0, ai > 0; (b). 



a2 > 0, ai < 0; (c). a2 < 0, ai > 0; (d). a2 < 0, ai < 0. 



Figure 7.4: Special orbits asymptotic to fixed points in the resonant annulus 



in forward or backward time when 



> 1, (a), ai > 0; (b). ai < 0. 
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Figure 7.5: Homoclinic and heteroclinic orbits for the case {N = 3, resonant 
region). 

7.4 Homoclinic and Heteroclinic Orbits for the Case {N = 3, 
Resonant Region) 

Inside the resonant annulus A, the dynamics for both the full system (7.2) 
and the first order system (7.3) is shown in figures 7.1 and 7.2. As a con- 
sequence of Theorem 5.1, there exist orbits asymptotic to the saddles in 
forward time. For fixed {e,a,uj), these orbits correspond to the intersection 
points between the surface (5.6) and the figure eight level sets of the saddles. 
In backward time, these orbits approach the resonant annulus A. Since the 
Hamiltonian is conserved and the figure eights are given as the level sets of 
the Hamiltonian, in backward time, these orbits also approach the saddles 
on the same figure eight level sets. Thereby, we have homoclinic or hete- 
roclinic orbits as shown in figure 7.5 for example. In drawing these figures, 
the fact that the quantity 



only changes 0(e) in finite time interval, and the fibers are smooth with 
respect to their base points and £. 



N-l 




n=0 
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8 Conclusion 



In this paper, we have proved the existence of homochnic tubes for the 
discrete cubic nonhnear Schrodinger equation under Hamiltonian perturba- 
tions, and discuss the symbohc dynamics of invariant tubes in the neigh- 
borhoods of such homochnic tubes when the system is four dimensional. 
In future works, we wih study such symbohc dyanmics structures when 
the dimension is higher. Since the system that we study here is a ncar- 
integrable Hamiltonian system, we are also interested in studying KAM tori 
and stochastic layer structures inside each invariant tube, in particular, the 
homochnic tube, and the connection with the symbolic dynamics of invari- 
ant tubes. Due to the normally hyperbolic nature of the center manifold, 
the relevant work on such KAM tori is that of Graff [5]. 

Acknowlegment: I am greatly indebted to Brenda Prazier for compil- 
ing all the artistic works. 
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